#ifndef RATE_LAW
#define RATE_LAW

double ReactionModel(double Xs, double P, double T);
double InvReactionModel(double w, double P, double T);
#endif